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1.  INTRODUCTION 
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The  research  described  In  this  report,  (cf.  previous  Interim  reports  [1, 
2,3]),  Is  generally  concerned  with  the  analysis  of  stress  In  the  vicinity  of 
cooling  holes  In  turbojet  engine  blades  with  transpiration  or  film  cooling. 

A schematic  diagram  of  a typical  blade1  is  shown  In  Figure  1.  Portions  of 
the  cross-sections  of  two  typical  types  of  blades  (taken  midway  between  cool- 
ing holes)  are  shown  In  Figures  2 and  3.  Possible  cooling  hole  arrangements 
In  the  vicinity  of  the  leading  edge  of  the  blade  are  as  Indicated  in  the  men- 
tioned Figures.  When  such  blades  are  in  service  the  blade  material  Is  sub- 
jected to  a thermomechanlcal  environment.  The  outer  blade  surfaces  are  In 
2 

contact  with  combustion  gases,  the  inner  surfaces  are  in  contact  with  cooling 
air,  and  spin  of  the  turbine  subjects  the  blade  to  a centrifugal  field.  Of 
particular  Interest  In  the  described  research  Is  the  "vicinity"  of  a single 
cooling  hole  (e.g.  Figure  4)  where  thermomechanical  stress  concentrations  of 
present  Interest  occur. 

The  method  of  analysis  employed  for  investigation  of  this  problem  Is  the 
Boundary  Integral  Equation  (BIE)  method  (cf.  [4,5]).  Under  assumptions  of 
linear,  homogeneous,  Isotropic,  steady-state,  thermoelastic  material  behavior, 
analytical  BIE  relations  applicable  to  an  arbitrary  portion  of  a blade  are 
written.  Next,  discretization  assumptions  are  made  leading  to  a discretized 
version  of  the  analytical  relations  In  the  form  of  simultaneous  linear 
algebraic  equations.  The  unknowns  In  these  equations,  generally  speaking,  are 
the  desired  quantities,  e.g.  stresses,  at  points  of  Interest.  Thus  a mathe- 


Speclflc  blade  configurations  and  data  supplied  by  General  Electric  (G.E.), 
Evendale,  Ohio. 
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All  data  In  this  report,  unless  otherwise  stated,  refer  to  first  stage  tur 
bine  blades  subjected  to  the  highest  temperatures. 
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matlcal  model  of  the  problem  Is  made  under  the  Indicated  constitutive  and 
discretization  assumptions.  Our  various  computer  programs  are  then  used  to 
provide  a numerical  solution  to  this  model  problem. 

Actually,  in  the  course  of  the  research  a number  of  specific  boundary 
value  problems  for  specific  facsimiles  of  portions  of  blades  are  posed  and 
solved.  Details,  rationale,  discussion,  and  specific  results  for  these 
problems  form  the  body  of  this  report  together  with  a discussion  of  the  latest 
version  of  the  BIE  analysis  capability  created  largely  for  and  used  In  our 
analyses. 

The  basic  integrity  of  our  working  computer  programs,  i.e.  the  reliability 
and  accuracy  to  be  expected  within  the  framework  of  the  theory  versus  the 
level  of  discretization  needed  for  the  numerical  analysis,  has  largely  been 
established  In  [1,2,3].  In  those  references  a number  of  representative  "test" 
problems  were  posed  and  solved  for  just  this  purpose.  In  the  present  report, 
however,  emphasis  is  placed  on  modelling  specific  problems,  the  results  from 
which  are  intended  to  give  Insight  into  mechanisms  for,  causes  of,  and  levels 
of  thermoelastic  stress  in  an  advanced  design  turbine  blade.  Accuracy  of 
results  as  a function  of  level  of  numerical  analysis  effort  Is  not  so  much 
the  Issue  here  as  Is  the  reality  and  reliability  of  the  models  themselves, 
especially  as  regards  synmetry  assumptions,  boundary  conditions,  and  thermal 
data  input  from  the  blade  manufacturer.  Discussion  of  these  matters  and  others 
is  included  in  this  report. 


2.  THE  BOUNDARY  INTEGRAL  EQUATION  METHOD 
FOR  THERMOELASTIC  PROBLEMS 


2.1  BRIEF  BACKGROUND 

The  BIE  Is  becoming  well  established  as  a companion  to,  or  In  some  cases, 

•j  a viable  alternative  for  more  conventional  methods  of  stress  analysis,  prin- 

cipal and  most  familiar  among  which  Is  the  Finite  Element  (FE)  method.  An 
extensive  background  for  the  BIE  or  a lengthy  a-priori  justification  for  its 
employment,  advantages,  etc.  In  problems  of  the  present  type  Is  out  of  place 
here.  Rather  we  call  attention  to  references  [4,9,15,16]  as  fairly  recent 
Indicators  of  some  of  the  features  of  the  BIE,  some  aspects  of  its  analytical 
foundations,  and  its  scope.  References  [1,2,3,5,13,14]  record  to  a large 
extent  the  progress  and  development  of  the  BIE  capability  as  used  in  the 
present  context  including  the  solution  of  a number  of  test  problems  some  of 
which  were  motivated  by  and  are  related  to  the  present  research. 

2.2  FORMULATION 

Consider  an  arbitrarily  shaped  solid  body  B in  a thermomechanical  environ- 
ment which,  as  a consequence  of  conditions  imposed  on  Its  surface  S and  poss- 

3 

Ibly  a body  force  distribution  b^  per  unit  volume  , suffers  a distribution  of 
temperature  0,  displacement  u^,  and  stress  from  a uniform  reference  state. 
The  stress,  the  displacement,  and  temperature  fields  throughout  the  body  are 
frequently  of  engineering  Interest. 

The  approach  to  such  problems  In  classical  thermoelastlclty  usually  In- 
volves two  steps:  (1)  Solve  a heat  transfer  problem  first  to  obtain  the 
steady  state  6 field  based  on  well-posed  thermal  boundary  conditions  on  S. 

I 

3 " ■■ 

Standard  Cartesian  Tensor  Analysis  Is  used.  Latin  subscripts  have  the 
range  1,2  or  1,2,3  according  as  B occupies,  respectively,  a region  of 
two  or  three-dimensional  space. 

L , . i m . _ 


- 
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(2)  Then  regard  6,^  as  a type  of  body  force  which,  together  with  well-posed 
mechanical  boundary  conditions,  contributes  to  the  u..  and  field  In  a 
familiar  way  (e.g.  [6]  p.  77),  and  solve  for  u..  and/or  a^. 

The  Boundary  Integral  Equation  (BIE)  method  for  such  thermoelastic 
problems  recognizes  the  above  two  steps  and  begins  with  an  integral  formula 
pertinent  to  the  heat  transfer  problem  in  the  form  (cf.  [4]  p.  1): 


c(p)6(p) 


[e'(q)K(p,Q)  - 0 (Q) K* (p,Q)]dS(Q). 


(1) 


In  formula  (1)  the  prime  (')  Indicates  the  (outward)  normal  derivative  of  the 
function,  K Is  a two-point-dependent  kernel  function  of  differing  form  (see 

O 

Appendix)  according  as  B pertains  to  a three-dimensional  (E  ) or  two-dimen- 
2 

sional  (E  ) body,  Q is  a point  e S,  and  p is  a point  e B or  S.  If  p e B, 

c(p)  * 4ir  or  c(p)  * 2it,  respectively,  according  as  the  body  B e E or  B e E . 

If  p c S,  c(p)  is  a number  dependent  only  upon  the  local  geometry  of  S at  p. 

In  practice  c(p),  p e S,  is  calculated  as  the  integral  of  K'  over  S.  This  is 

apparent  from  setting  6 = 1 in  formula  (1).  If  S has  a unique  tangent  at  p, 

3 2 

c(p)  ■ 2tt  for  B e E and  c(p)  * it  for  B t E as  expected. 

The  significance  of  formula  (1)  is  as  follows.  Let  p e S and  regard  (1) 
as  a relation  among  6 and  6'  on  S pertaining  to  one  and  the  same  temperature 
field  throughout  B.  Now  since  only  part  of  the  pair  (0,0')  (e.g. 0 or  0'  but 
not  both)  may  be  specified  at  any  Q e S according  to  a well-posed,  steady 
state,  heat  transfer  problem,  formula  (1)  with  p e S provides.  In  principle, 
a means  of  obtaining  that  part  of  the  pair  (0,0')  not  Initially  prescribed. 

That  is,  we  regard  formula  (1)  with  p e S as  an  Integral  equation  to  be  solved 
for  the  unprescribed  parts  of  the  pair  (0,0')  on  S.  Once  (0,0')  are  completely 
known  everywhere  on  S,  consider  p e B In  (1),  and  that  formula  yields  0 at  any 
point  In  the  Interior  of  B In  terms  of  the  Indicated  Integral  over  S of  known 


Liu 


quantities.  Thus,  assuming  the  above-described  process  can  be  carried  out, 
and  the  method  for  doing  so  numerically  will  be  Indicated  subsequently. 

Step  (1)  In  classical  thermoelasticity  has  been  completed. 

Step  (2)  Is  similar  to  Step  (1)  and  involves  the  following  formula 
(cf.  [5])  pertinent  to  the  thermoelastic  problem: 

,(Q)^(Q)n1(Q)}U1j(p,Q)-u<(Q)T^(p,Q)]dS(Q) 


c1j(p)ui(p)  » 


U WU  WIC  l 

fs  [{V 


‘1’ 


ij' 


+ Js  C{tKQ)+Y9(Q)}W,,;j(p,Q)  - {^(Q)+Ye(Q)},W,j(p,Q)]dS(Q) 


+ k. 


L- 


(Q)W(p,Q)dS(Q). 


(2) 


In  formula  (2)  U^,  and  W are  two-point-dependent  kernel  functions  of 
differing  forms  (see  Appendix)  according  as  B e E or  6 e E , t^(Q)  Is  the 
surface  traction,  r(p,Q)  = |Q-p|,  (Q)  are  unit  (outward)  normal  components, 

4»(Q)  Is  the  surface  value  of  a function  i|>  defined  throughout  B such  that  the 
mechanical  body  force,  if  any,  is  derivable  from  ij>  according  to 


bi  * -i| »ti  with  ^,i1  = k£ 


(3) 


where  kQ  Is  a constant,  as  Is  y (see  Appendix).  In  formula  (2)  p may  be  c B 


or  S and  the  (')  signifies  normal  derivative  as  before.  If  p e S,  c^(p)  Is 


a matrix  of  numbers  determined  In  a manner  similar  to  that  for  determining 
c(p)  In  formula  (1)  (see  e.g.  [5],  Eqs.  (1),  (14),  (15)).  If  p e B, 


c1j*p)  s 61j* 


The  significance  of  formula  (2)  Is  as  follows.  Let  p e S and  regard  (2) 
as  a relation  among  u^  and  t^  on  S and  the  functions  6,  e',  i|i,  <|>'  also  on  S, 
all  of  which  pertain  to  one  and  the  same  thermoelastic  stress  state  throughout 
B.  Now  6 and  61  are  known  on  S following  completion  of  step  1 as  previously 
described.  (Note  that  setting  p e B and  solving  for  Interior  0 with  formula 
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(1)  is  not  required  in  using  formula  (2)).  Variables  i|>  and  t|>‘  are  regarded 
as  known  boundary  quantities  obtainable  from  a potential  for  real  mechanical 
body  force.  What  is  not  known  completely,  of  course,  is  the  pair  (u^.t^)  on 
S,  since  only  part  of  the  pair  can  be  prescribed  corresponding  to  well-posed 
boundary  conditions.  For  example,  u^  may  be  prescribed  on  S (displacement 
type  problem),  or  t^  may  be  prescribed  on  S (traction  type  problem),  or  parts 
of  each  may  be  prescribed  over  portions  of  S (various  kinds  of  mixed  problems), 
etc.  Thus,  the  solution  strategy  via  formula  (2)  with  p e S is  to  solve  the 
integral  equation  for  that  part  of  the  pair  (u^,t^)  not  initially  prescribed. 
Once  (u.,t^)  are  completely  known  everywhere  on  S,  consider  p e B in  (2)  and 
that  formula  yields  u^  at  any  point  In  the  interior  of  B in  terms  of  the 
indicated  integrals  over  S of  known  quantities.  Finally,  the  stress  components 
o^  at  any  p e B are  given  via  formula  (2)  by  first  obtaining  the  displacement 
gradient  u^,k  at  p e B.  This  is  done  simply  by  differentiating  the  kernel 
function  expressions  with  respect  to  x^  at  p and  evaluating  the  indicated  in- 
tegrals of  the  differentiated  kernel  functions.  Based  on  Hooke's  Law,  the 
stress  field  is  given  In  terms  of  u^.j  according  to 

a1j  3 T^T  Vm6ij  + vCu1’j  + uj*i]  ’ 61jY0  (4) 


wherein  p e B.  If  on  S is  desired,  it  may  be  obtained  directly  in  terms 


of  (6»Uj,t^)  on  S (actually  e.t^  and  "in-surface"  derivative(s)  of  u^  on  S are 


required)  in  a familiar  fashion  (see  e.g.  [7,8]).  Thus  for  surface  o^.  there 


Is  no  need  to  refer  to  the  interior  p. 

This  completes  the  analytical  formulation  of  the  two-step  thermoelastic 
solution  process  via  the  BIE  method. 

2.3  NUMERICAL  PROCEDURES 


To  Implement  the  solution  strategy  outlined  above,  we  first  discretize  the 
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surface  S of  B,  B e E , into  a total  number  m of  surface  elements  of  either 
quadrilateral  or  triangular  shape  as  shown  in  Figure  5(a)  and  locate  eight  or 
six  nodes,  respectively,  on  each  element.  The  cartesian  coordinates  x“  of 
each  node  are  specified,  and  the  cartesian  coordinates  of  a non-nodal  point 
of  an  element  are  assumed  to  be  given  by  (cf.  [5,9]) 


x,(S)  ■ 


a = 1 ,2. . . ,6  or  8 


(5) 


in  which  m“(e)  are  second-order  shape  functions  of  intrinsic  coordinates 
(0  = We  use  one  set  ^or  quadrilateral  shapes  and  another  set 

for  triangular  shapes  (see  Appendix).  Next,  we  assume  that  any  of  the  functions 
of  Q alone  in  equations  (1)  and  (2)  varies  over  the  elements  according  to 

4>(C)  - M0^)*01 


(6) 


where  <pa  is  the  nodal  value  of  the  particular  variable  at  issue  and  the  shape 
functions  are  identical  to  those  used  in  equation  (5).  Based  on  the  discre- 
tizing assumptions,  formula  (1)  becomes 


= a.  eta 


c(PTi)e(PTi)  = 0 


I 


- 0' 


ota 


L 


M(5)K(PT1,Q(5))J(C)d5 


M(OK'(Pn,Q(C))J(OdC 


(7) 


and  formula  (2)  becomes 


‘ijC 


n»“l<pn>  " nf  Js 

■ ““  Is 


Mct(5)Uij(Pn,Q(C))J(5)d5 
Sa 


Ma(0T1j(Pri,Q(5))J(C)dC  + T™  I ^(€)W,,j(Pn,Q(0)J(C)d5 


- T' 


ota 


I, 


Ma(5)W,i(Pn,Q(0)J(0dC  + k? 
J n U(J*1 


L ■ 


(8) 


j(Q(0)W(Pri,Q(C))J(0dC 


f1 
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In  which  It  is  assumed  p = P e S,  and  in  which  n.,  * T = \p  + y0; 


and  where  the  quantities  ^ have  the  meaning  "the  value  of  $ at  local  node  o 


on  surface  element  S ".  Note  that  a - 1,2,...,  m and  n * 1,2,...,  n,  where 


n Is  the  total  number  of  nodes;  summation  Is  implied  on  all  repeated  sub- 
scripts; J(C)  is  simply  the  ratio  d$(Q)/d5,  which  more  explicitly  is 

J(0  - |a.,  x a2|  (9) 


where  the  components  of  the  vectors  in  the  cross  product  are  given  by 


aui 


- 3^*  ..a 


35  Xi* 
y 


(10) 


Note  that  since  a node  always  lies  on  at  least  two  surface  elements  and 
frequently  more,  is  the  same  quantity  for  several  sets  of  values  of  a and 
o.  This  little  bookkeeping  problem,  when  addressed,  leads  to  viewing  the  sys- 
tem of  equations  (7)  and  (8)  in  the  following  matrix  forms,  respectively: 

F6  = 66'  * (11) 


Au  - Bfl  * CT  - OT'  + E. 


(12) 


In  equation  (11)  the  column  matrices  6 and  6'  are  n elements  long  and 
contain  the  elements  0®,  0,e,  8 * 1,2,...  n,  which  are,  respectively,  temp- 
erature and  normal  derivative  of  temperature  at  global  node  number  8.  The 
square  matrices  F and  G are  n x n in  size  and  contain  the  Integral  coefficients 


(and  the  c(P  ) term)  as  Indicated  in  equation  (7). 

n 


In  equation  (12),  the  column  matrices  u and  J1  are  3n  elements  long  and 

fi  ^ 


contain  the  elements  u!j,  8 a 1,2,...,  n,  which  are,  respectively,  displace- 
ment and  "augmented"  traction  components  at  global  number  8;  T and  T'  are  n 


fi  fi 

elements  long  and  contain  T , T , i.e.,  the  surface  contribution  of  temperature 


and  body  force  potential  and  normal  derivatives  of  these  quantities,  as  in- 
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dicated  above,  at  the  nodes.  The  square  matrices  A and  B are  3n  x 3n  elements 

mm  *v 

In  size  and  contain  the  Integral  coefficients  of  u and  £1  as  Indicated  In 

**m  mm 

equation  (8).  Similar  remarks  hold  for  the  n x 3n  size  matrices  C and  D. 

mm  mm 

The  column  E contains  the  3n  values  of  the  last  Integral  In  equation  (8)  as 
Indicated. 

The  previously  outlined  solution  strategy  when  referred  to  equations  (11) 
and  (12)  above  amounts  to  the  following:  Discretize  the  surface  S of  B Into  a 
number  m of  surface  elements  and  form  and  store  the  matrices  F and  G.  Specify 
an  appropriate  combination  of  elements  6 and  0'  at  the  nodes  corresponding  to 
a well-posed  heat  transfer  problem.  Perform  the  indicated  matrix  multipli- 
cations in  (11)  involving  specified  quantities  and  solve  the  system  (11)  for 
the  elements  of  6,  6'  not  previously  specified.  At  this  point,  the  boundary 

•M  rnm 

4 

set  0,  8*  are  completely  known. 

mm  mm 

To  begin  Step  (2),  form  and  store  the  matrices  A,  B,  C,  D,  E.  Specify 

mm  mm  mm  mm  mm 

the  elements  T and  T'  which  Involve  the  now-known  nodal  values  of  the  set 

mm  mm 

0,  0'  and  also  (assumed)  known  nodal  values  of  the  body  force  potential  i|> 

mm  mm 

and  its  normal  derivative  ip' . Next,  perform  the  indicated  matrix  multipli- 
cations to  yield  a single  known  vector  on  the  right  side  of  equation  (12). 

Then  specify  a sufficient  number  of  elements  u and  t corresponding  to  a well- 

mm  mm 

posed  elastostatlcs  problem.  Perform  the  remaining  Indicated  matrix  multi- 
plications In  (12)  Involving  known  or  specified  quantities  at  this  stage,  and 
solve  the  system  (12)  for  the  remaining  unknowns.  Now  the  set  u,  t is  com- 

mm  mm 

pletely  known  such  that  the  full  complement  of  surface  quantities  u,  t,  0,  0', 

mm  mm  mm  mm 

\p,  xp*  pertaining  to  the  desired  thermoelastic  stress  state  In  B is  available 

mm  mm 

at  the  nodes. 

I 

If  desired,  although  not  necessary  for  step  two  In  the  analysis,  we 
could  take  p c B In  formula  (7)  and  generate  the  temperature  0 at  any 
desired  Interior  p. 


.... — 


^ 


Putting  p e B f In  formula  (8)  now  would  yield  u^(p)  at  any  desired  p. 
ai j (?)  would  be  obtained  similarly  by  first  obtaining  ^ (p)  from  formula  (8) 
as  previously  indicated  and  then  using  Hooke's  Law  (4).  Boundary  are  ob- 
tainable In  terms  of  0,  u,  t as  indicated  earlier. 

** 

2 

For  B e E , i.e.  a plane  problem,  S Is  a curve  or  the  union  of  several 
curves,  and  most  of  the  preceding  remarks  apply  intact  with  the  appropriate 
reduction  in  dimension.  Note  that  Sa  is  now  an  element  of  boundary  curve  with 
nodes  at  each  end  and  one  (usually  midway)  between.  There  are  now  three  shape 
functions  M“U),  a * 1,2,3,  which  have  the  form  shown  in  the  Appendix,  and  £ 
is  a single  scalar  variable.  Matrices  which  were  formerly  3n  or  3n  x 3n  in 
size  are  now  2n  or  2n  x 2n,  respectively. 

Having  outlined  the  solution  process  for  thermoelastic  boundary  value 
problems  via  the  BIE  method  above  in  some  detail,  it  is  apparent  that  success 
with  the  method  depends.  In  large  part,  on  being  able  to  obtain  sufficiently 
accurate  values  for  the  elements  of  the  matrices  A,  B,  C,  D,  E,  F,  G.  Exam- 
ination of  equations  (7)  and  (8)  reveals  that  the  mentioned  matrix  elements 


are  sums  of  Integrals  of  the  form 


[“  5 f 

v k 


rf‘(C)r(l>n,Q(C))0(C)dC 


In  which  Sa  Is  a square  or  equilateral  triangle  as  shown  In  Figure  5(b)  (or 
the  interval  -1  to  1 of  the  straight  line  for  B e E ),  since  a mapping  to  those 
shapes  Is  Inherent  In  the  change  of  variables  Implied  by  the  introduction  of 
the  shape  functions  (5).  r Is  any  of  the  tensor,  vector,  or  scalar  kernel 
functions  Indicated  In  equations  (7)  and  (8).  The  quantities  c(P  ),  c^ i(P  ) 

In  equations  (7)  and  (8)  respectively  prove  also  to  be  the  sum  of  certain 
Integrals  of  the  type  (13)  (see  [5],  Eqs.  (14)  and  (15)  for  details).  Thus, 
it  is  necessary  to  obtain  a sufficiently  accurate  value  for  all  l“a  by  some 
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numerical  Integration  scheme.  For  a given  choice  of  surface  discretization 
(number  m and  arrangement  of  surface  elements  and  number  of  nodes  n),  the  man- 
ner of  obtaining  values  for  Is  perhaps  the  most  important  issue  In  the 
solution  process.  Total  computer  (CPU)  time  for  a given  problem  Is  largely 
a function  of  how  this  Is  done. 

3 

To  expand  on  this  matter  for  B e E , note  first  that  two  distinct  cases 
In  the  evaluation  of  l“a  arise:  (a)  P^  t Sa,  and  (b)  P^  e Sa<  In  case  (b), 
note  that  since  certain  of  the  kernel  functions  are  0(l/r  ) and  others  at 
least  0(l/r),  the  Integrands,  while  Integrable,  require  special  attention. 

The  first  case,  (a).  Is  more  straightforward.  For  Integration  over  the 
square  we  use  a Gaussian  quadrature  formula  with  weight  function  1.0  to  obtain 

l“  » / F°(P  ,5)d£  4 I1  S2  A(V  A1 Fa(Pn,E(Xl)5(X2))  (14) 

na  JSa  n p*l  V=1  w v 11 

where 

Fa(Pn.O  i Mfx(E)r(Pri,QU))J(5)  (15) 

and  A ^l\  A^*2^  are  weight  factors  corresponding  to  the  order  Xj,  X2  chosen 
for  the  quadrature  formula  [10].  For  Integration  over  the  triangle  we  use  a 
special  triangularly  symmetric  quadrature  rule  by  Lyness  [11],  which  symbol- 
ically may  be  written  In  the  same  form  as  (14).  In  either  case,  accuracy  Is 
In  principle  a function  of  the  order  (size  of  Xj,  X2)  of  the  quadrature  rule 
but  so  also  Is  the  cost  (CPU  time),  so  that  proper  choice  of  the  pair  X^,  X2 
Is  an  important  practical  Issue  for  a given  Integration.  After  considerable 
study  and  experimentation  (details  will  be  available  In  [12])  we  have  auto- 
mated the  choice  of  X^,  X2  as  a function  primarily  of  the  ratio  of  the  dis- 
tance of  the  point  Pn  to  some  characteristic  size  of  the  element  Sff.  Second- 
arily, we  Increase  the  order  X^  X2  for  a relatively  "distorted"  element  over 
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and  above  the  order  otherwise  called  for  by  the  "distance"  criterion.  Both 
criteria  are  an  attempt  to  choose  Xj,  X2  just  large  enough  for  a given  inte- 
grand to  give  values  of  Integrals,  via  formula  (14),  with  consistent  error. 

For  case  (b),  where  P^  e Sa,  we  first  take  account  of  the  fact  that 

c(P  ) and  c.JP  ) are  calculated  In  a manner  (cf.  [5],  Eqs.  (14)  and  (15)) 
n ■ J n 

2 

which  analytically  requires,  for  r of  0(l/r  ),  that  we  evaluate 

L 06) 

* o 

with  6 * 1 whenever  node  a = P ; otherwise,  6 * 0.  Thus,  if  we  imagine  a 

an  n an 

local  polar  coordinate  system  p,4>  as  shown  In  Figure  (6)  for  typical  locations 
of  P^,  we  note  that  tf*(5)  Is  0(p)  for  a ^ P^  and  that  (M^O-l)  is  of  the  same 
order  for  a * P^.  Thus  since  rfP^.E)  is  0(l/p  ) and  d£  becomes  pdpd<f> • the  In- 
tegrand in  (16)  is  no  longer  Improper  and  may  be  evaluated  using  the  same 
types  of  quadrature  rules  previously  Indicated. 

In  case  (b),  since  P„  occupies  one  of  the  nodes,  the  "distance"  criteria 

fl 

for  choice  of  X-| , X2  no  longer  has  meaning  and  a threshold  choice  is  made 

according  to  a preselected  accuracy  requirement.  Finally,  as  with  the 

P^  i Sa  case,  we  can  increase  the  order  X^,  X2  for  a relatively  "distorted" 

element  if  necessary. 

2 

For  B e E , for  P^  t So,  a one-dimensional  version  of  formula  (14)  is 

used.  Here,  the  choice  of  order  of  Integration  (size  of  X)  as  a function  of 

the  proximity  of  P to  S Is  less  critical  than  In  three  dimensions  because  of 
no 

the  comparatively  fewer  total  Integrations  to  be  done.  Hence,  we  presently 
pick  a threshold  size  for  X (probably  too  high  for  most  Integrations)  and 
compute  all  P^  t Sc  Integrals  for  constant  X.  Efficiency  measures,  neverthe- 
less, could  be  Introduced  as  for  three  dimensions.  Regarding  the  P^  e 
case,  note  that  rf*(0  - 6^  Is  0(5)  and  since  the  kernel  functions  are  at 


4 


I 
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most  0(1/5)  no  transformation  comparable  to  the  polar  coordinate  transformation 
Is  necessary  in  two  dimensions.  There  Is  the  problem,  however,  of  one  inte- 
grand of  0(log  5)  which  can  be  integrated  easily  enough  using  a "log"  weighted 
Gauss  or  by  removing  the  singularity  (cf.  [7],  [13],  [14]). 

Based  on  the  foregoing  discussion,  numerical  values  of  all  elements  of 
the  coefficient  matrices  A,  B,  etc.  can  be  calculated  automatically  by  a 
digital  computer  upon  specifying  the  cartesian  coordinates  of  all  the  nodes 
on  the  surface  S of  a body  under  investigation,  the  material  parameters,  and 
the  quadrature  parameters  X-j,  Xg. 

All  subsequent  computations  were  performed  on  an  IBM  370-165  digital 
computer  using  single-precision  arithmetic.  The  program  and  all  data  for  a 
given  problem  (Step  (1)  or  (2))  were  retained  in  core.  Note  that  only  the 
matrix  of  coefficients  of  unknowns  at  the  nodes  need  be  formed  and  stored  in 
a square  array.  The  other  elements  of  matrices  of  coefficients  may  be 
Immediately  multiplied  by  other  quantities  as  they  are  generated  and  stored 
In  column  vectors.  Individual  coefficient  arrays  need  not  be  formed  and 
stored  separately  as  Implied  earlier. 

2.4  FOUR  COMPUTER  PROGRAMS 

Based  on  the  formulation  In  Section  2.2  using  the  discretization/numerical 
procedures  of  Section  2.3,  four  separate  computer  programs  have  been  created. 
Two  of  the  programs  perform  heat  transfer  analyses  (Step  1)  for  problems  In 
two  and  three  spatial  dimensions,  respectively.  The  remaining  two  programs 
similarly  perform  thermoelastic  analyses  (Step  2).  These  four  programs  form 
the  working  tools  for  the  research  to  be  described,  and  they  have  been  devel- 
oped, checked  and  exercised  In  various  stages  (cf.  [1,2,3])  throughout  the 
course  of  AFOSR  sponsorship  of  the  present  project. 

|! 
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3.  THE  PRIMARY  RESEARCH  PROBLEM  DEFINED 
3.1  THE  BLADE  ENVIRONMENT 

When  in  service,  the  outer  surface  of  a typical  turbine  blade  (see 
Figure  1)  is  in  contact  with  hot  combustion  gases  with  nominal  temperatures 
of  the  order  2400°F.  As  indicated  in  Figures  2 and  3,  the  blade  has  hollow 
chambers  through  which  cooling  air,  at  a nominal  temperature  of  1400°F,  is 
pumped  under  pressure  and  spewed  out  through  the  cooling  ("gill")  holes  and 
vent  holes  (cf.  Figures  1 through  4).  One  effect  of  this  air  flow  through 
the  hollow  blade  interior  is  to  lower  the  average  blade  temperature.  Further, 
a cooling  film  of  air  is  provided  over  the  blade  surface  in  the  vicinity  of 
the  gill  holes.  This  is  of  particular  importance  at  the  leading  edge  of  the 
blade  where  temperatures  are  highest.  The  flow  of  air  and  hot  gases  and  the 
centrifugal  field  provide  the  thermomechanical  environment  for  the  blade 
material.  All  flows  and  fields  are  assumed  steady  in  the  present  analysis. 

Clearly,  from  the  above  description,  stress  is  induced  in  the  blade 
material  as  a consequence  of  (a)  the  nonuniform  temperature  distribution 
from  place  to  place  in  the  blade,  (b)  the  centrifugal  field  on  the  blade, 
and  (c)  the  gas  (and  air)  pressure  and  drag  acting  as  surface  loads  or 
tractions  on  the  blade.  Each  of  these  causes  of  stress  may  be  regarded 
separately  and  the  effects  superposed.  However,  according  to  engineering 
analysts  at  General  Electric,  cause  of  stress  (c)  leads  to  significant  values 
only  near  the  throat  (lower  portions)  of  the  blade  and  such  values  are  known 
to  be  only  a fraction  of  that  expected  due  to  cause  of  stress  (b).  Cause 
of  stress  (a),  the  nonuniformity  of  the  temperature  field  is  the  least  under- 
stood and  forms  the  thrust  of  the  present  research. 

The  thermal  stress  problem  is  especially  interesting  due  to  stress 
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raising  character  of  the  gill  holes.  These  holes  contribute  to  elevated 
stress  in  the  vicinity  potentially  in  two  ways:  (1)  the  nominal  thermal 
stress  field  due  to  the  gross  blade  temperature  distribution  is  concentrated 
near  the  hole  surfaces  due  to  the  very  geometrical  presence  of  the  holes, 
and  (11)  the  holes  provide  for  localized  additional  variation  of  temperature 
due  to  cooling  effect  of  the  air  carried  by  the  holes.  As  will  be  shown, 
the  geometry  effect  (1)  is  far  more  important  than  (11).  Therefore  for 
purposes  of  the  present  research  we  consider  first  a fixed  blade  subjected 
to  zero  traction  and  body  force  but  subjected  to  a nonuniform  temperature 
field  to  be  described  in  some  detail.  Later  we  superpose,  in  appropriate 
fashion,  the  effect  due  to  the  steady  centrifugal  field. 

3.2  THE  PRIMARY  RESEARCH  PROBLEM  MODEL 

We  are  most  interested  in  the  stress  in  the  immediate  vicinity  of  a 
single,  typical  gill  hole.  Our  first  model  for  the  "vicinity"  of  such  a 
hole  for  a blade  with  a cross-section  of  the  type  in  Figure  2 is  as  shown 
In  Figure  4.  This  Is  a thick  irregular  platelike  body  obtained  by  "cutting" 
along  the  dotted  lines  of  Figures  2 and  1.  This  body  is  bounded  as  indi- 
cated In  Figure  4 by  portions  of  the  circular  "vertical"  cooling  chamber 
surfaces,  the  three  "vertical"  cut  surfaces,  the  "front  vertical"  surface, 
the  "top  horizontal"  cut  surface  located  midway  between  two  typical  gill 
holes,  the  "bottom  horizontal"  cut  surface  through  a typical  gill  hole,  and 
finally  by  the  "half"  gill  hole  surface  Itself.  The  angle  of  inclination  a 
of  the  gill  hole  with  the  front  surface  is  variable  but  for  the  specific 
body  of  Figure  4 and  our  subsequent  analysis  of  that  particular  body. 


a - 60°. 

Of  special  Interest  now  in  the  body  of  Figure  4 is  the  distribution  of 


normal  stress  along  the  edge  AB  of  the  gill  hole.  Of  particular  Interest  is 
the  peak  stress  which  occurs,  as  might  be  expected,  at  the  "feather  edge"  of 
the  hole,  point  A. 

To  solve  for  the  stress(es)  of  interest  via  the  analysis  method  outlined 
j in  Section  2,  it  is  necessary  to  first  solve  the  heat  conduction  problem  for 

the  boundary  of  the  body.  This  requires  a proper  statement  of  well -posed 
thermal  boundary  conditions  on  each  surface  of  the  body  of  Figure  4.  Follow- 
ing completion  of  this  first  step,  9 and  0'  (cf.  Section  2.2)  will  be  known 
everywhere  on  all  surfaces  S.  The  second  step  which  uses  as  input  the  now 
known  pair  (0,0')  on  S,  requires  a proper  statement  of  well -posed  mechanical 
boundary  conditions  on  each  of  the  surfaces.  Output  from  the  second  step 
will  include  stresses  along  the  edge  AB  and  in  particular  at  point  A. 

Clearly,  the  analysis  as  outlined  for  the  body  of  Figure  4 depends 
heavily  on  a proper  statement  of  the  thermal  and  mechanical  boundary  condi- 
tions on  each  surface.  Hence  the  rationale  for  these  boundary  conditions 
which  in  turn  leads  to  the  consideration  of  other  problem  models  will  now  be 
considered  in  some  detail. 

3.3  BOUNDARY  CONDITIONS 

The  blade  portion  of  the  body  in  Figure  1 is  not  quite  prismatic  but  it 
is  nearly  so  and  the  assumption  of  identically  shaped  cross-sections  taken 
normal  to  the  "vertical"  or  z direction  over  a blade  with  Indefinite  length 
(z  direction)  is  convenient  and  easily  justifiable  for  present  purposes. 
Moreover,  thermal  conditions  are  assumed  to  be  periodic  functions  of  z with 
a period  of  one  gill  hole  spacing  such  that  the  fields  of  Interest,  e.g. 

c 

stress,  displacement,  are  also  similarly  periodic  . The  described  periodicity 
r 

The  centrifugal  field  Is  assumed  zero  et.  seq.  until  stated  other- 
wise. 
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or  "cyclic  symmetry"  clearly  demands  the  following:  (a)  the  heat  flux  0' 
is  zero  on  the  horizontal  cut  surfaces  of  Figure  4,  (b)  each  of  these  sur- 
faces must  remain  plane  under  thermal  load,  and  (c)  the  shear  traction  must 
be  zero  on  each  of  these  surfaces. 

Now  from  typical  thermal  data  supplied  by  General  Electric,  it  is  pos- 
sible to  supply  a non-zero  value  of  0‘  over  the  vertical  cut  surfaces  of 
Figure  4.  Unfortunately,  the  variation  of  such  0'  with  z,  while  probably 
not  great  over  the  period,  is  unknown  and  we  thus  assume,  in  the  absence  of 
better  information,  0'  to  be  independent  of  z on  the  mentioned  surfaces. 

For  the  remaining  surfaces,  i.e.  the  vertical  part-cylindrical  chambers, 
the  gill  half-hole  surface,  and  the  front,  gas/gas-air  mix,  contacting  sur- 
face, typical  data  for  k,  hc  and  0Q  in  a convective  boundary  condition  of 
the  type 

(k/hc)o*  * (0Q-0)  (17) 

were  made  available  by  General  Electric.  In  Equation  (17),  k is  the  thermal- 

conductivity,  h Is  a convective  surface  heat  transfer  coefficient  and  0 is 
c 0 

a "contact-fluid"  temperature.  Numerical  values  for  h and  0 over  the 

v V/ 

relevant  surfaces  were  made  available  such  that  with  the  mentioned  known 
values  of  01  on  the  cut  surfaces.  It  Is  possible  to  pose  and  solve  a heat 
transfer  problem  for  the  body  of  Figure  4.  Again,  hc  and  0Q  are  taken  In- 
dependent of  z.  Nevertheless,  the  presence  of  the  gill  holes  yields  0 and 
0*  which  depend  on  z over  the  gas  or  air  contacting  surfaces  as  required. 

With  the  thermal  boundary  data  in  hand  it  remains  to  consider  the 
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mechanical  boundary  conditions  on  the  surfaces  of  Figure  4.  Now  as  men- 
tioned earlier  we  may  assume  the  gas  or  air  contacting  surfaces  to  be 
traction  free.  Thus  the  gill  half-hole  surface,  the  vertical  part-cylind 
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rlcal  surfaces  and  the  front  surface  have  the  simple  condition  t = 0. 

On  the  lower  cut  horizontal  surface  we  may  state  the  symnetry  condition  as 
t$  * 0 and  un  * 0 where  s and  n refer  to  in  plane  and  normal  directions, 
respectively.  On  the  upper  cut  horizontal  surface,  the  condition  that  it 
remain  plane  and  shear  traction  free  is  helpful  but  not  as  immediately 
complete  as  the  lower  surface  condition  since  the  location  of  the  plane  re- 
mains unknown.  We  will  return  to  this  boundary  condition  shortly. 

The  vertical  cut  surfaces  of  Figure  4 present  the  most  interesting 
consideration  of  boundary  conditions.  The  tractions  on  these  surfaces 
clearly  arise  from  the  surrounding  portions  of  the  blade  from  which  the 
piece  is  cut.  These  tractions  are  non-zero  to  the  extent  that  the  thermal 

B 

field  in  the  whole  blade  is  non-linear  in  the  cartesian  coordinates.  Thus 
to  obtain  representative  values  for  the  tractions  t on  these  vertical  cut 
surfaces  it  is  necessary  to  consider  gross  blade  behavior.  We  accomplish 
this  in  the  following  way. 

The  thermal  conditions  as  already  noted  are  periodic  in  the  variable  z, 
and  under  the  prismatic  assumptions,  this  periodic  variability  is  strictly 
a consequence  of  the  presence  of  the  gill  holes.  In  their  absence  the  ther- 
mal field  would  be  plane  or  independent  of  z.  The  question  arises  then  to 
what  extent  do  the  gill  holes  alter  the  thermal  tractions  which  would  exist 
on  the  vertical  cut  surfaces  of  Figure  4 in  the  absence  of  gill  holes  ass- 
uming, of  course,  comparable  overall  temperature  distributions  in  the  blade 
in  the  two  cases?  This  question  is  sensibly  addressed  as  follows.  Examine 
the  cross-section  of  Figure  2.  Ignore  the  geometrical  influence  of  the 
periodically  arranged  gill  holes  and  perform  a heat  transfer  (Step  1)  and  a 
subsequent  thermoelastic  analysis  (Step  2)  on  the  cross-section  of  the  blade 


shown  bounded  externally  by  the  airfoil  shape  and  internally  by  the  circular 
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boundaries.  For  this  plane  (strain)  analysis  use  actual  thermal  boundary 
conditions  for  that  portion  of  the  blade  in  Figure  2 shown  in  Figure  7. 
Typical  numerical  data  for  the  boundary  conditions  are  available  from  the 
manufacturer.  Use  also  the  (in-plane)  zero  traction  condition  and  calculate 
the  tractions  on  the  three  dotted  lines  of  Figure  2,  which  are  the  edges  of 
the  three  vertical  cut  surfaces  of  Figure  4.  These  computed  thermal 
tractions  are  clearly  an  approximation,  nevertheless  an  excellent  approxi- 
mation we  think,  to  the  actual  tractions  existing  on  the  vertical  cut 
surfaces  of  Figure  4.  The  difference  between  these  two  sets  of  tractions 
would  be  due  to  the  influence  of  the  cooling  holes  on  the  gross  thermo- 
mechanical behavior  of  the  blade  which  gave  rise  to  the  computed  tractions 
in  the  manner  just  indicated.  We  argue  that  the  difference  would  be  slight 
especially  on  the  plane  midway  between  gill  holes  and  the  tractions  computed 
as  indicated  represent  an  excellent  approximation  to  the  "true"  tractions. 

The  latter  would  be  clearly  unavailable  without  an  incredibly  complicated 
three  dimensional  analysis  of  a body  with  the  same  cross-section  as  in 
Figure  2 but  at  least  several  gill  holes  thick  in  which  the  holes,  despite 
their  small  size  and  probable  influence  on  the  variables  at  present  issue, 
would  have  to  be  modelled. 

On  the  assumption  that  the  tractions  computed  via  the  above  described 
planar  analysis  represent  a valid  set  to  be  placed  on  the  vertical  cut 
surfaces  of  Figure  4,  it  now  remains  to  precisely  state  a boundary  condition 
for  the  top  cut  horizontal  surface  argued  earlier  to  remain  plane  and  shear 
traction  free. 

Toward  this  end  note  first  that  the  magnitude  of  normal  traction 
tn  = tz  (t$  ■ 0)  acting  in  the  out-of-  plane  z direction,  under  the  plane 
strain  assumptions  for  the  plane  thermoelastic  solution  for  the  cross-section 
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of  Figure  7,  is  given  by 


*z  ■ v(oxx  + “yy*  • “E9 


(18) 


wherein  axx  and  o are  the  in-plane  normal  stresses  available  from  the  men- 


tioned plane  strain  solution.  Numerical  values  of  tz  are  thus  available. 


in  practice,  at  a selected  array  of  points  on  the  material  occupied  area  A 
of  Figure  7. 

We  may  now  imagine  a prismatic  body  of  indefinite  length  with  the  cross- 
section  of  Figure  7 (no  gill  holes)  under  the  plane  0(x,y)  distribution,  and 
traction  free  on  its  lateral  surfaces.  On  any  cross-section  (i.e.  Figure  7) 


t acts  and  the  deformation  in  the  body  is  plane.  If  we  now  calculate  the 


resultant  axial  force  F on  this  body  and  the  components  of  resultant  couple 


according  to 
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and  define  tractions  tQ  and  t^  respectively  given  by 
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(19) 
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we  may  superpose  to  yield  a resultant  traction. 


t mto*h  +lz 


(22) 


In  the  above,  xy  are  coordinates  relative  to  the  centroid  of  the  cross-section 

of  Figure  7,  and  I , I , I are  the  area  moments  and  product  of  inertia 
xx  yy  xy 

of  A relative  to  the  xy  coordinates.  Clearly,  the  net  force  and  couple 
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associated  with  t are  zero.  Further,  the  prismatic  body  subjected  to 
* 6 

0(x,y)  and  t only  is  a representative  model  of  a typical  length  of  the 
blade  in  a static  thermoelastic  state  ignoring  end  effects,  the  gill  holes, 
and,  of  course,  the  centrifugal  field. 

It  is  a relatively  straightforward  analytical  exercise  to  now  calculate 
the  displacement  of  one  cross-section  relative  to  another  due  to  the  t + t^ 
distribution  acting  on  this  bar.  In  particular  we  can  find  the  displacement 
of  one  cross-section  one-half  gill  hole  spacing  away  from  another.  Such  a 
cross-section,  like  all  others,  remains  plane  under  tQ  + t-j  load.  The  posi- 
tion of  this  plane,  as  well  as  others,  is  clearly  a function  of  the  gross 
temperature  field  in  the  blade  uninfluenced  to  any  significant  effect  by  the 
actual  presence  of  the  gill  holes  which  occupy  only  an  insignificant  portion 
of  the  volume  of  the  blade  material.  Note  that  conditions  of  symmetry  demand 
that  even  in  the  presence  of  gill  holes,  cross-sections  through  gill  holes 
must  remain  plane.  Clearly  these  planes,  with  or  without  gill  holes,  are  al- 
most identical.  Thus  we  now  have  our  boundary  condition  for  the  top  surface 
of  the  body  of  Figure  4,  namely,  a normal  displacement  un  f 0 which  is  identi- 
cal with  the  mentioned  analytically  determined  displacement  arising  from  the 
tractions  tQ  and  t, . 

Note  that  knowledge  of  the  non-zero  displacement  un,  computed  as  described, 

on  the  top  surface  of  Figure  4 is  a more  meaningful  boundary  condition  than  a 
* 

prescription  of  t . This  is  of  the  utmost  significance  in  the  present  analysis 

* 

and  is  true  for  the  following  reason.  While  the  t actually  present  over  the 
small  part  of  the  cross-section  (top  surface  of  Fig.  4)  may  be  significantly  diff- 
erent than  that  given  by  Equation  (22) , i.e.  locally  altered  by  the  presence  of  the  gill 

Note  that  Figure  7 depicts  an  incomplete  blade  cross-section.  The  "tail" 
portion  is  missing.  A brief  discussion  of  the  effect  of  this  on  the  analy- 
sis is  given  on  page  27  paragraph  2. 
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hole,  un  is  computed  from  tQ  and  which  in  turn  are  based  on  the  distribu 

tion  of  t_  over  the  entire  area  of  Figure  7.  Thus  the  u_  used  is  substan- 
z n 

tially  the  same  with  or  without  gill  holes. 


We  now  have  a complete  set  of  well -posed  mechanical  boundary  conditions 
for  the  thermoelastic  body  of  Figure  4.  Numerical  results  and  further  dis- 
cussion of  this  "Primary  Research  Problem  Model"  are  given  in  Section  4. 

3.4  THE  EFFECT  OF  THE  CENTRIFUGAL  FIELD 


By  comparison  to  the  above  consideration  of  boundary  conditions,  the 
boundary  conditions  on  the  body  of  Figure  4 to  take  account  of  the  cen- 
trifugal field  of  the  blade  are  fairly  simple.  Using  the  condition 
yn  = 0,  t$  = 0 on  the  lower  horizontal  cut  surface,  tp  = oQ,  t$  = 0 on  the 
upper  horizontal  cut  surface,  with  all  other  surfaces  traction  free,  is 
reasonable  and  consistent  with  the  manufacturer's  practice.  Isothermal 
conditions  are  assumed  and  the  presence  of  centrifugal  "body"  force  over 
the  small  volume  itself  may  be  sensibly  neglected  by  comparison  to  aQ  which 
approximates  the  centrifugal  pull  of  the  adjacent  portion  of  the  blade  across 
the  surface  on  which  oQ  acts.  Again  the  distribution  of  normal  stress  along 
the  edge  AB  of  the  gill  hole  is  the  desired  output  from  the  solution  of  this 
problem. 

Note  that  the  condition  u = 0 on  the  lower  horizontal  surface  of 

n 

Figure  4 assumes  that  this  surface  is  a plane  of  symnetry  for  the  body  plus 
its  reflection  in  that  surface.  This  body  according  to  the  preceding  state- 
ments is  loaded  in  simple  uniform  tension  on  its  upper  and  lower  surfaces. 


How  well  this  approximates  the  actual  state  of  stress  due  to  the  centrifugal 
field  is  somewhat  a function  of  the  gill  hole  spacing  vs.  gill  hole  diameter. 
Indeed,  if  the  gill  holes  are  sufficiently  far  apart  the  assumption  of  uni- 
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form  a0  midway  between  holes  is  a good  one.  However,  in  the  extreme,  where 
the  holes  are  very  close  together,  e.g.  less  than  one  hole  diameter  apart, 
the  above  assumption  is  probably  poor.  In  the  latter  case  a better  approxi- 
mation would  be  to  first  recognize  that  the  surfaces  between  holes  while  not 
plane  are  similarly  deformed  under  centrifugal  loading  and  it  is  perhaps  not 
too  gross  an  assumption  to  take  the  portions  of  the  surfaces  involved,  i.e. 
the  top  surface  of  Figure  4,  its  mirror  image,  and  therefore  also  the  bottom 
horizontal  surface  of  Figure  4 to  be  parallel  planes.  Thus,  an  alternative 
boundary  condition  for  the  body  of  Figure  4 to  simulate  centrifugal  loading 

would  be  u = 0,  t = 0 on  the  lower  surface,  u = u , t = 0 on  the  upper 
•i  v n o «> 

surface. 

Numerical  data  for  unit  input  values  of  oQ  andu^,  thereby  using  both 
assumptions,  are  given  and  compared  in  Section  4. 
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4.  SOLVED  PROBLEMS  FOR  THE  PRIMARY  MODEL 
4.1  THE  PLANE  PROBLEMS  FOR  FIGURE  2 

According  to  the  discussion  of  boundary  conditions  in  Section  3.3  it  is 
necessary  to  first  solve  the  heat  conduction  problem  for  the  blade  cross-sec- 
tion of  Figure  7.  Values  of  hc  and  0Q  were  given  by  G.E.  on  a scale  drawing 
from  which  we  were  able  to  state  a thermal  boundary  condition  of  the  type  (17) 
for  each  part  of  the  boundary  shown  in  Figure  7.  Note  that  the  gill  hole, 
being  below  or  above  the  cross-section  of  Figure  7,  is  not  part  of  the  bound- 
ary for  this  analysis  which  is  planar.  Using  the  discretization  depicted  in 
Figure  7 with  m = 67  line  segments  and  n = 134  nodes,  our  heat  conduction 
analysis  yielded,  for  the  specific  boundary  data  used,  the  isotherms  shown  in 
the  figure.  The  depicted  temperature  distribution  compares  favorably  with  the 
manufacturer's  temperatures  obtained  by  other  means. 

Next,  assuming  in-plane  tractions  on  the  body  of  Figure  7 to  be  zero,  we 
calculate,  through  Step  (2)  of  the  analysis,  the  stress  distributions  across 
the  three  surfaces,  i.e.  the  three  dotted  lines  in  Figure  7.  Results  and 
material  properties  are  shown  in  Figure  8.  Distribution  of  normal  stress  is 
shown  In  detail.  Directions  of  shear  stresses,  which  are  an  order  of  magni- 
tude  smaller  than  the  normal  stresses,  are  indicated  by  a single  arrow  as 
shown.  Thus,  the  major  effect  of  the  thermal  field  in  the  blade  cross-section 
is  to  produce  bending,  via  the  three  "thermal  couples"  shown  in  Figure  8 on 
the  piece  of  the  blade  indicated.  While  we  have  only  sketchy  stress  data  from 
the  manufacturer  with  which  to  compare,  our  results  are  believable  when  com- 
pared  with  the  essential  features  of  thermal  stress  in  a hollow  cylinder,  cool 
on  the  interior  and  hot  on  the  exterior.  This  problem  was  treated  in  some 
detail  In  [2,5].  Additional  discussion  of  the  cylinder  problem,  and  how  a 
major  stress  inducing  consequence  of  the  temperature  distribution  in  the 


blade  material  is  present  in  the  cylinder  problem  will  be  given  subsequently. 

The  distribution  of  stress  shown  in  Figure  8,  assumed  to  be  uniform  over 

the  length  of  the  vertical  cut  surfaces  of  Figure  2,  represents  the  needed 

tractions  on  the  mentioned  surfaces,  as  discussed  in  Section  3.3.  These 

tractions  represent  the  major  mechanical  influence  of  the  surrounding  blade 

material  on  the  "piece"  in  Figure  2 as  required.  The  remaining  influence, 

i.e.  the  traction  t is  readily  obtained.  A plot  of  the  t distribution  is 

shown  in  Figure  9.  Approximately  245  points  were  positioned  over  the  area  A 

of  Figure  2 at  which  a , a , and  e were  obtained  from  which  t was  computed 

xx  yy  z 

via  Equation  (18).  From  this  t distribution,  F,  C , and  C from  Equation  (19) 

z ^ y 

were  obtained.  Based  on  a nominal  constant  reference  temperature  of  1800°F, 
the  computed  t2  distribution  yielded  the  following  values 

F = -498.60  lb,  C = 3.03  in-lb,  C * -45.19  in-lb  (23) 

x y 

Therefore  tQ  = 9,082  lb/in2  from  Equation  (20).  Further,  A = .055  in2,  Ixx  * 

2.77  x 10"4  in4,  I = 7.96  x 10’4  in4,  I = 0.11  x 10"4  in4.  From  these 
yy  xy 

numerical  data  the  linear  t^  distribution  is  now  known  via  Equation  (21).  Thus 

the  plane  displacement  un  = unQ  + u^  of  any  cross-section  of  the  type  in 

Figure  7 relative  to  any  other  is  now  obtainable  as  the  simple  deformation  due 

to  the  tractions  t + t|.  In  particular  the  displacement  un  of  the  "top" 

surface  of  Figure  4 relative  to  the  "bottom"  surface  is  available  and  could 

be  used,  as  indicated  in  Section  3.3  as  our  "top"  surface  boundary  condition. 

However,  top  surface  boundary  condition  may  be  simplified  further  by 

noting  that  the  actual  values  of  F,  C , and  C listed  above  are  somewhat 

x y 

arbitrary.  Specifically,  it  is  clear  that  any  change  in  the  arbitrarily 
chosen  reference  temperature  of  1800°F  would  lead  to  a proportional  change  in 
tQ  or  F.  However,  the  corresponding  constant  part  of  the  displacement  un. 
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L., 


i.e.  unQ,  would  be  similarly  changed  and  the  arbitrariness  compensated  in  the 


superposition  process.  Similar  remarks  apply  for  t. , C , C , and  the  linear 

i x y 

part  of  the  displacement  uR,  i.e.  unj,  since,  in  principle,  an  arbitrary  non- 
constant but  linearly  varying  reference  temperature  distribution  could  have 
been  chosen.  Indeed,  there  exists  a Capri ori  unknown)  linear  reference 


temperature  distribution  for  which  F,  C , and  C would  all  be  zero.  For  this 

x y 


(non-arbitrary)  reference  distribution,  the  condition  ur  = 0 on  the  "top" 


surface  of  Figure  4 would  be  the  proper  boundary  condition.  This  linear 
reference  distribution  may  now  be  obtained  as  follows.  Note  first  that 


u = o £ = 1.12  x 10"5  in. 
no  t 


(24) 


wherein  tQ  = 9,082  lb/in2,  E = 24.3  x 106  lb/in2,  and  l = 0.03  in  is  one-half 


the  gill  hole  spacing  in  the  present  design  configuration.  Secondly,  since 


I is  an  order  of  magnitude  smaller  than  I and  I and  since  C is 

xy  aa  yy  a 


similarly  small  compared  to  Cy,  the  bending  implied  by  Equation  (21)  is  nearly 


symmetrical  (i.e.  the  neutral  surface  of  bending  is  inclined  only  8.8°  from 
the  y axis  (Figure  7)).  Thus  with  negligible  error  for  present  purposes  t^ 
is  approximately  given  by 

'yy 


C x 

ti 4 s ”5*68  x 1q4  x 1b/in2 


wherein  C. 


-45.19  lb-in  and  I = 7.96  x 10"4  in4. 


form 
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-7.01  x 10  J x in. 


(25) 

Therefore,  u^  has  the 

(26) 
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Now  the  desired  reference  temperature  distribution  has  the  form 


0R  s a + bx  + cy 


(27) 


J 
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where  a,  b,  and  c are  constants.  Since  the  unknown  u has  the  form 

n 

un  = a(a  + bx  + cy)l  » (1.12  - 7.01x)  x 10-5  in.  (28) 

in  which  a * 9.2  x 10"®/°F  is  the  coefficient  of  thermal  expansion,  it  is 
clear  that  the  choices 

a = 40.6°F  + 1800°F,  b = -254.0°F/in,  c = 0 (29) 

gives  the  explicit  form  for  the  reference  distribution.  Subsequent  analyses 

for  the  primary  research  problem  model  of  Figure  4,  based  on  this  particular 

eR  as  the  reference  temperature  distribution  may  use  the  simple  top  surface 

boundary  condition  * o.  All  arbitrariness  with  regard  to  reference 

temperatures  in  subsequent  calculations  has  thereby  been  removed. 

There  is  one  remaining  difficulty  with  the  modelling  and  superposition 

solution  process  under  the  present  circumstances.  Unfortunately,  no  thermal 

or  other  data  compatible  with  the  data  for  the  portion  of  the  blade  cross- 

section  shown  in  Figure  7 were  available  for  the  "tail"  portion  not  shown. 

To  have  made  a more  representative  model,  tQ  and  t-j  should  be  computed  on  the 

basis  of  F,  C , and  C (obtained  from  t ) for  the  enti re  cross-section.  How- 
* y z 

ever,  modification  of  our  present  analysis,  should  data  for  the  tail  portion 
of  the  blade  become  available,  is  readily  made.  A new,  perhaps  more  repre- 
sentative, value  for  un  could  easily  be  computed  analytically  leading  to  a 
(possibly)  new  reference  temperature  distribution.  We  suspect  alterations  in 
the  following  results  would  be  slight. 

Having  obtained  all  the  necessary  data  from  the  above  described  plane 
analysis,  we  proceed  now  to  the  three-dimensional  analysis  for  the  primary 
model  (cf.  Figure  4,  discussion  in  Section  3.3). 

4.2  HEAT  CONDUCTION  ANALYSIS 

Figure  10  depicts  the  primary  research  problem  model  (cf.  Figure  4)  with 
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dimensions  to  scale  and  a partial  view  of  the  discretization  pattern  (ref. 

Sec.  2.3)  used  for  this  part  of  the  analysis.  The  number  of  surface  elements 
m * 69  and  the  number  of  nodes  n * 215.  The  circular  cooling  chambers  have 
diameters  of  60.0  mils  and  140.0  mils,  respectively,  the  gill  hole  diameter 
is  16.5  mils,  and  the  length  l (z  direction)  of  the  piece  is  30.0  mils,  i.e., 
one-half  the  gill  hole  spacing  for  this  particular  configuration.  These 
dimensions  and  discretization  pattern  are  used  for  all  subsequent  thermoelas- 
tic analyses  of  the  primary  research  problem  model. 

As  discussed  in  the  second  paragraph  of  Section  3.3  we  now  input  the  data 
supplied  by  G.E.  to  Equation  (17)  for  each  of  the  fluid  contacting  surfaces. 
Part  of  these  data  were  applicable  to  the  previously  performed  planar  analysis 
for  the  cross-section  of  Figure  7,  but  now  we  have  the  gill  hole  and  the  length 
l to  contend  with  and  thus,  of  course,  a three-dimensional  problem.  Despite 
the  presently  necessary  assumptions  about  0'  not  varying  with  z,  and  similar 
assumptions  on  hc  and  0Q  in  Equation  (17)  (ref.  parag.  2 Sect.  3.3),  a heat 
conduction  analysis  for  the  body  of  Figure  10  under  the  stated  boundary 
conditions  leads  to  a smooth  and  believable  three-dimensional  variation  of 
temperature  over  the  body.  Since  the  subsequently  relevant  outcome  of  this 
heat  conduction  analysis,  i.e.  0 and  01  at  all  the  nodes  on  the  surface  of 
the  body  are  not  of  intrinsic  Interest,  and  difficult  to  display  meaningfully, 
these  data  are  not  shown.  No  data  from  General  Electric  with  which  to  com- 
pare any  of  our  three-dimensional  temperature  variations  were  available,  of 
course,  which  is  consistent  with  z varying  input  data  being  unavailable  from 
G.E.  as  previously  mentioned.  Nevertheless,  our  solution  in  its  x,y  variation 
is  consistent  with  their  data  and  our  own  plane  solution  previously  described. 
Further,  the  local  thermal  effect  of  the  gill  hole  manifests  itself  realisti- 
cally by  cooling  the  piece  of  Figure  10  by  amounts  and  in  places  that  are  con- 
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sistent  with  reasonable  expectations.  Step  (1)  for  the  research  problem  model 
is  now  completed. 

4.3  THE  FULL  THERMOELASTIC  PROBLEM 

With  the  complete  set  (6,9')  known  numerically  at  all  the  nodes  for  the 
body  of  Figure  10,  we  begin  Step  (2)  by  supplying  the  mechanical  boundary 
conditions  as  follows.  The  tractions  of  Figure  8,  assumed  constant  (with 
respect  to  z)  over  the  "vertical"  cut  surfaces  of  Figure  10,  are  input  as 
numerical  values  of  normal  and  shear  traction  at  appropriate  nodes.  Values 
of  traction  at  nodes  on  all  fluid  contacting  surfaces  are  zero.  Following 
the  adjustment  in  (0,0')  to  provide  the  specific  reference  temperature  dis- 
tribution 0R  (ref.  Eqs.  (27)  - (29)),  we  use  the  condition  un  = 0,  t$  = 0 
at  all  nodes  on  the  "top"  horizontal  surface.  This  condition  is  now 
identical  with  that  on  the  "bottom"  surface  of  Figure  1C. 

Our  thermoelastic  solution  based  on  the  stated  boundary  condition  yields 
the  thermoelastic  stress  field  throughout  the  body  of  Figure  10.  Of  primary 
interest  is  the  normal  stress  distribution  on  the  bottom  surface,  (cf.  Figure 
4)  along  the  line  AB,  i.e.  one  edge  of  the  hole.  Point  A is  of  special  in- 
terest, i.e.  the  "feather  edge".  The  stress  distribution  along  this  line  and 
along  the  "front  edge"  AC  are  shown  in  Figure  11.  Further,  the  "hoop  stress" 
distribution,  i.e.  the  stress  at  the  surface  of  the  hole  in  directions  tan- 
gent to  the  hole,  is  shown  in  Figure  12. 

The  stress  at  point  A as  shown  is  tensile,  and  our  analysis  predicts,  as 
probably  expected,  a high  gradient  of  stress  as  point  A is  approached  along 
lines  AB  and  AC.  Compression  is  indicated  along  most  of  the  length  AC  whereas 
tension  is  indicated  along  the  entire  length  AB7.  Higher  levels  of  tensile 

7 These  results  disagree  with  data  reported  for  a related  problem  in  [3], 
Figure  7.  An  error  was  discovered  which  invalidates  that  data  in  [3]. 
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stress  at  the  inside  (point  B)  end  of  the  gill  hole  are  indicated  than  at  the 
feather  edge  (point  A).  Further,  the  magnitude  of  compressive  stress  at  C is 
comparable  to  that  of  the  tensile  stress  at  A.  These  results  require  some 
interpretation  and  discussion. 

As  an  aid  to  interpretation  of  our  numerical  data,  to  help  assess  their 
reliability,  and  to  gain  greater  insight  into  the  mechanism  for  the  inducement 
of  stress  in  the  body  of  Figure  10,  we  imagine  the  following  separate  problems 
which,  when  "added  together",  yield  the  stress  results  depicted  in  Figures  11 
and  12. 

Consider  the  body  of  Figure  10  at  uniform  temperature  but  loaded  on  its 
"vertical"  cut  surfaces  with  the  tractions  of  Figure  8.  The  top  is  traction 
free  (unconstrained),  and  the  bottc;n  (et.  seq.)  is  held  plane  by  normal  trac- 
tions only,  which,  as  before,  are  the  quantities  of  interest.  This  problem, 
as  far  as  the  body  of  interest  is  concerned,  is  isothermal,  but  the  load  is 
a set  of  tractions  which  arise  from  the  fact  that  the  body  is  part  of  and, 
of  course,  acted  upon  by  the  surrounding  blade  - the  whole  being  in  a thermal 
environment.  Stress  distributions  along  lines  AB  and  AC  for  this  problem 
(called  component  problem  (i),  et.  seq.)  are  shown  in  Figure  13,  and  the  hoop 
stress  distribution  is  shown  in  Figure  14. 

The  peak  tensile  stress  at  point  A for  this  loading  is  about  twice  the 
value  shown  in  Figure  11.  The  way  the  tractions  are  applied  (cf.  Figure  8) 
clearly  is  consistent  with  high  tensile  stress  at  A and,  less  intuitively 
perhaps,  consistent  with  the  hoop  stress  pattern.  This  means,  of  course,  that 
the  other  "components"  of  this  problem,  namely,  (ii)  the  actual  6/0  distribu- 
tion over  the  body  of  Figure  10,  and/or  (iii)  the  constraint  condition  un  = 0 
at  the  top  surface  provide  sufficient  compression  to  reduce  the  tensile  stress 
at  point  A by  50%  and  otherwise  yield  the  dominant  compression  along  most  of 

j 


31 


1 


line  AC,  as  shown  in  Figure  11.  Further,  these  other  "components"  must  be 
responsible  for  the  high  tension  along  the  line  AB. 

To  assess  the  separate  roles  of  components  (ii)  and  (iii),  we  consider 
the  body  of  Figure  10  with  no  applied  tractions  whatsoever  except  for  the 
bottom  which  has  the  same  boundary  condition  as  before.  Here  we  put  in  the 
’!  0^0  distribution  (referred  to  0R)  and  thus  solve  the  component  problem  (ii) 

to  assess  the  influence  of  the  (non-linear)  temperature  distribution  on  the 
piece  of  Figure  10,  removing  the  constraining  influence  of  the  surrounding 
material.  Here,  we  find  that  the  levels  of  stress  along  lines  AB,  AC,  and 
the  hoop  stress  distribution  are  comparatively  insignificant.  Our  solution 
of  this  problem  yields  nominal  (compressive)  stresses  of  the  order  two  to 
three  thousand  pounds  per  square  inch  alono  line  AC,  for  example,  with  an 
increase  to  about  ten  thousand  right  at  point  A.  The  stress  level  along  line 
AB  is  similarly  low.  Thus  it  seems  clear  that  the  constraint  or  mechanical 

effect  (albeit  arising  from  thermal  conditions  in  the  overall  blade)  of  the 

' 

blade  material  on  the  piece  of  Figure  10  dominates  in  producing  the  stresses 

< 

of  interest.  That  is,  component  problems  (i)  and  (iii)  dominate.  More 
specifically,  the  lateral  tractions  of  Figure  8 lead  to  high  tensile  stress  at 
the  feather  edge  (point  A),  which  is  reduced  by  the  compression  induced  by  the 
constraint  (un  = 0)  on  the  top  surface  of  Figure  10,  and,  to  a lesser  extent, 
by  the  0 ^ 0 distribution  on  the  body  itself.  Further  discussion  of  this  and 
related  matters,  and  the  assumptions  on  which  the  various  described  problems 
are  based  will  be  made  subsequently. 

4.4  STRESS  DUE  TO  THE  CENTRIFUGAL  FIELD 

As  mentioned  in  Section  3.4,  we  consider  two  models,  i.e.  two  problems, 
each  with  a different  top  surface  (Figure  10)  boundary  condition,  to  represent 
the  centrifugal  field. 

I I ^ 
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The  first  model  assumes  t = o = 1 on  the  top  surface  of  Figure  10  to 

no 

simulate  the  centrifugal  field.  With  this  model  and  the  previously  used 
discretization  pattern,  our  isothermal  elastic  analysis  yields  the  stress 
distributions  along  the  edge  of  the  gill  hole  AB  and  the  front  edge  AC  as 
shown  in  Figure  15.  The  maximum  tensile  stress  as  shown  in  Figure  15  yields 
a stress  concentration  ^max/°0  = 9.4.  This  figure  is  in  satisfactory  agree- 
ment with  the  stress  concentration  of  8 for  the  related  problem  of  a flat 
plate  with  a 60°  inclined  hole  in  tension  (discussed  at  length  in  [2]), 
considering,  of  course,  the  more  irregular  shape  of  our  body. 

The  second  model  assumes  up  = uQ  = 1 on  the  top  surface  of  Figure  10  to 
simulate  the  centrifugal  field.  As  noted  in  Section  3.4,  this  boundary  con- 
dition is  probably  better  than  the  uniform  aQ  for  closely  spaced  gill  holes. 
The  present  spacing  of  60  mils  is  typical,  although  we  understand  there  exist 
configurations  of  from  considerably  less  than  60  mils  spacing  up  to  120  mils 
or  greater  for  a gill  hole  diameter  of  16.5  mils.  Really,  any  spacing/dia- 
meter ratio  could  be  expected,  and  the  two  types  of  boundary  conditions  used 
above  probably  provide  upper  and  lower  bounds  on  the  centrifugally  induced 
stress-concentration  for  the  presently  used  ratio  of  hole  spacing  to  diameter. 
The  same  is  likely  true  for  any  other  ratio. 

Once  again  our  isothermal  elastic  analysis  yields  the  "edge  AB"  and 
"edge  AC"  distributions  as  shown  in  Figure  16.  Here  the  displacement  uQ  = 1 
is  statically  equivalent  to  an  average  applied  tensile  traction  a . Based 

o ” y 

on  oaVg,  a stress  concentration  at  the  feather  edge  of  the  hole  (point  A, 
Figure  16)  is  a /a  „ = 5.7.  This  figure  is  considerably  lower  than  the 
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9.4  above,  as  is  to  be  expected  from  the  inherent  difference  in  the  type  of 
input  (cf.  discussion  [2]  pg.  19). 
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5.  DISCUSSION 

5.1  ADEQUACY  OF  DISCRETIZATION 

The  question  of  integrity  of  the  numerical  data  presented  in  the  previous 
sections  is  largely  a function  (a)  of  the  adequacy  of  the  discretization  used 
for  the  various  problem  models  and  (b)  the  significance  and  reliability  of  the 
boundary  conditions,  symmetry  assumptions,  etc.,  and  input  thermal  data.  Both 
of  these  issues  are  subject  to  some  reexamination  as  follows. 

The  discretizations  used  for  the  plane  problems  are  more  than  adequate  to 
the  task  at  hand.  Our  experience  with  such  plane  discretizations,  and  the  out- 
come of  our  experiments  in  this  regard  with  the  present  problems,  leaves  little 
doubt  that  the  results  of  our  plane  analyses  are  as  accurate  as  the  knowledge 
of  input  thermal  data  and  stated  boundary  conditions  would  justify.  The  dis- 
cretizations used  for  the  three-dimensional  problems  admittedly  could  be  finer, 
to  provide  more  resolution  in  the  regions  of  highest  stress  gradients,  for 
example.  However,  finer  three-dimensional  discretizations  come  only  at  consid- 
erable cost  and  effort  (cf.  discussions  in  [1,2,3]).  Moreover,  the  work 
[1,2,3]  was  heavily  devoted  to  assessing  the  reliability  of  discretizations 
comparably  fine  to  those  used  here  in  connection  with  problems  for  which  we 
could  judge  accuracy  directly.  Indeed,  those  problems  were  chosen  for  their 
essential  similarity  to  the  ones  at  hand.  Further,  those  problems  helped  to 
exercise  and  check  all  of  the  important  thermal  and  mechanical  loading  effects 
and  geometrical  features  of  the  solved  problems  in  Section  4.  For  example, 
the  various  isothermal  stress-concentration  problems  solved  in  [1,2],  with 
comparably  fine  or  cruder  discretizations,  lend  confidence  to  that  part  of  the 
stress  of  interest  due  to  mechanical  loading  on  the  "top"  surface  of  Figure  10. 
Confidence  in  the  temperature  distributions  and  thermally  induced  stresses  is 
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similarly  gained  from  the  work  reported  in  [2],  One  '‘loading  feature"  in  the 
present  analysis,  not  tested  in  [1,2],  is  that  of  the  body  of  Figure  10  loaded 
essentially  in  bending  by  the  "couples"  of  Figure  8 leading  to  the  results  of 
Figures  13  and  14.  Thus,  to  lend  the  level  of  confidence  in  this  part  of  the 
analysis  that  we  have  in  the  other  parts  as  a consequence  of  earlier  work,  we 
consider  the  problem  described  in  the  next  section. 

5.2  PLATE  WITH  STRAIGHT-THRU-HOLE  IN  BENDING 

Figure  17  shows  to  scale  one  quarter  of  an  elastic  plate  with  dimensions 
10  x 16  x 4.  The  radius  of  the  hole  is  unity,  and  the  discretization  pattern 
used  to  solve  the  problem  to  be  posed  is  as  indicated  in  the  figure.  We  load 
this  plate  in  simple  lateral  bending  by  applying  a linearly  distributed  bending 
traction,  the  tensile  part  of  which  is  as  shown  in  Figure  18.  The  output  for 
this  problem,  the  stress  distribution  (tensile  part  only)  along  the  mid-surface 
through  the  hole, is  shown  in  the  same  figure.  The  stress  concentration,  de- 
fined as  the  ratio  of  the  peak  stress  at  the  hole  to  the  nominal  bending  stress 
from  our  solution,  is  2.49.  This  compares  favorably  with  an  expected  value 
from  the  literature  of  2.62.  Thus  with  the  discretization  pattern  shown,  which 
is  comparably  fine  to  that  shown  in  Figure  10  for  the  research  problem,  we  are 
able  to  achieve  peak  stress  accuracy  to  within  5%. 

Of  course,  the  hole  of  Figure  10  is  sharply  inclined  to  the  "front"  sur- 
face, which  leads  to  higher  stress  gradients  and  places  more  demand  on  the  dis- 
cretization pattern  than  the  straight-thru  hole.  Nevertheless,  judging  from 
the  predicted  stress  concentrations  near  straight-thru  holes  under  bending 
and  tension,  the  bending  stress  concentration  is  of  the  same  order  of  magnitude 
as  the  tensile  one.  Thus  since  we  have  already  (cf.  [2])  tested  a comparably 
fine  discretization  for  the  inclined  hole  under  tensile  load  we  have  reason 
to  believe  that  peak  stress  data  in  Figures  11  through  14  are  as  accurate  as 


L 


35 


I 


-r 

the  data  in  Figures  15  and  16.  Even  an  accuracy  estimate  of  perhaps  10%,  or 
a more  pessimistic  estimate  of  20%  to  25%  accuracy,  in  our  judgement,  would 
still  render  the  mentioned  data  of  value  in  providing  insight  into  the  general 
picture  of  how  stress  is  induced  in  the  vicinity  of  the  gill  holes. 


5.3  SIGNIFICANCE  OF  BOUNDARY  CONDITIONS 


The  question  raised  earlier  of  the  significance  and  reliability  of  the 
assumptions  of  cyclic  symmetry  on  the  turbine  blade,  assumed  to  be  prismatic 
and  indefinitely  long,  and  the  necessary  consequences  of  plane  surfaces  should 
now  be  reexamined  in  light  of  the  results  for  component  problems  (i),  (ii), 
and  (iii)  in  Section  4.3. 

Recall  that  the  problems  are  (i)  the  "isothermal"  body  under  the  "lateral 
couples"  of  Figure  8,  (ii)  the  "free"  body  subjected  to  only  the  9^0  dis- 
tribution, and  the  implied  third  problem  (iii),  which  is  the  body  of  Figure  10 
subjected  to  the  top  surface  condition  un  = 0.  In  the  three  problems,  we  sub- 
mit that  the  boundary  condition  un  = 0 on  the  top  surface  for  problem  (iii)  is 
most  questionable.  This  in  turn  questions  the  validity  of  the  un  = 0 condi- 
tion on  the  full  thermoelastic  problem  as  originally  posed  and,  therefore,  the 
validity  of  the  data  in  Figures  11  and  12. 

If  the  blade  were  very  long  (we  assumed  indefinitely  long),  the  stated 
boundary  condition  and  thus  the  numerical  data  (Figures  11  and  12)  would  be 
reasonably  representative  for  gill  holes  away  from  the  ends.  On  the  other 
hand,  for  a gill-hole  near  the  end  of  the  blade,  the  top  surface  boundary 
condition  more  properly  should  be  t = 0 as  was  used  in  component  problem  (i) 
(Figures  13  and  14).  Thus,  ignoring  the  contribution  of  component  problem 
(ii),  it  appears  that  the  tensile  stress  at  the  feather  edge,  point  A,  lies 
somewhere  between  the  two  extremes.  That  is  the  43,000  figure  is  probably 
too  high  for  most  gill  holes,  and  the  19,400  figure  is  probably  too  low. 
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especially  for  the  holes  nearer  the  free  end.  Of  course,  the  lower  gill  holes 
are  subject  to  the  larger  stress-concentration  due  to  the  centrifugal  field, 
so  there  is  a compensating  effect,  i.e.  higher  tensile  stress  due  to  the  ther- 
mal field  where  the  centrifugal  effect  is  smaller  and  vice-versa. 

It  is  possible  that  if  a finer  picture  of  the  thermal  distribution  over 
the  surfaces  of  Figure  10  were  available  from  the  manufacturer,  especially  as 
regards  the  variation  with  depth  of  the  parameters  hQ  and  0Q,  component  pro- 
blem (ii)  could  have  made  a more  significant  contribution  to  the  stress  fields 
of  interest.  However,  our  basic  finding  that  the  stresses  of  component  pro- 
blem (i),  i.e.  the  stresses  induced  by  the  lateral  traction  distribution  :• 
most  important  is  consistent  with  the  information  available  from  an  idealized 
but  essentially  similar  problem. 

Tim  idealized  problem  consists  of  a thick  circular  cylinder  with  cool 
inner  surface  and  hotter  outer  surface  held  in  a configuration  of  plane  strain. 

This  problem  was  discussed  at  length  in  [2]  pg.  9 ff,  and,  upon  reflection, 
it  exhibits  many  of  the  important  characteristics  of  the  primary  research 
problem.  We  have  two  cool  inner  surfaces  (circular  cooling  chambers)  instead 
of  one,  but  the  outer  surface  is  hot  in  both  cases.  The  transverse  stress 
in  both  cases  (cf.  [2],  Table  I and  the  present  Figure  8)  produces  bending 
on  a portion  of  the  cylinder  as  it  does  on  our  research  body.  Moreover,  if 
one  considers  only  a piece  of  the  circular  cylinder,  e.g.  a piece  subtending 
an  angle  of  say  m/4  radians,  and  subjected  to  the  lateral  bending  tractions 
from  the  rest  of  the  cylinder  contacting  the  piece,  the  inplane  stresses  in- 
duced by  those  tractions  under  isothermal  conditions  for  the  tt/4  piece  are 
indistinguishable  from  those  obtained  if  the  actual  temperature  distribution 
over  the  piece  is  included.  This  observation  lends  confidence  to  our  conclu- 
sion that  component  problem  (ii)  is  relatively  insignificant.  In  short,  the 
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transverse  stress  of  Figure  8 acting  on  the  body  of  Figure  10  is  essentially 
similar  to  that  in  the  cylinder  problem.  This  is  the  stress  that  "is  concen- 
trated by"  the  gill  hole  and  of  primary  interest.  This  observation  is  contrary 
to  our  earlier  expectations  [1,2,3]  that  perhaps  the  local  9 field  on  the 
piece  of  Figure  10  was  at  least  equally  important.  Thus,  the  problem  of  a 
circular  cylinder,  cool  inside  and  hot  outside,  with  a hole  drilled  from  outer 
to  inner  surface  (at  a variable  angle  to  the  vertical  direction)  would  be  an 
excellent  model  for  the  important  phenomena  in  the  gill-hole-cooled  turbine 
blade.  This  model  is  especially  applicable,  we  judge,  for  the  blade  cross- 
section  of  Figure  3.  If  this  is  so,  the  bending  problem  of  the  previous 
section  has  additional  significance  in  the  following  sense.  Imagine  the 
applied  couples  to  be  thermally  induced,  since  the  plate  is  imagined  cut  from 
the  blade  of  Figu*^  3.  Then  the  isothermal  analysis,  previously  argued  to 
be  valid  for  the  small  piece,  yields  the  stress  concentration  of  about  2.5. 
Variations  on  this  problem,  i.e.  a parameter  study, could  be  performed, 
by  varying  the  angle  of  inclination  of  the  gill  holes  and,  although  less  im- 
portant perhaps,  allowing  the  plate  to  have  some  curvature.  These  problems 
would  clearly  simulate  the  conditions  in  the  vicinity  of  the  gill  holes  de- 
picted in  Figure  3.  Such  a parameter  study  could  be  readily  done  with  our 
developed  computing  capability  pending  the  necessary  thermal  data  for  specific 
blade  configurations  from  the  blade  manufacturer. 

A final  reference  to  the  cylinder  problem  is  in  order.  This  problem  was 
solved  in  [2]  under  plane-strain  assumptions.  If  we  imagine  a superposed  uni- 
form traction  tQ,  acting  on  the  ends  of  the  cylinder,  sufficient  to  render  the 
cylinder  axial-force-free  (or  equivalently  use  a proper  reference  temperature 
0^),  we  have  an  essentially  similar  situation  to  that  of  the  primary  research 
body  of  Figure  10.  Specifically,  a segment  of  the  cylinder  is  under  an  axial 
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traction  t which  is  tensile  near  the  inner  surface  of  the  cylinder  and  com- 
pressive near  the  outer  surface.  This  pattern  is  consistent  with  the  data  of 

★ 

Figure  11.  Further,  Figure  19  depicts  the  t distribution  on  the  top  surface 
of  the  research  body,  and  it  will  be  noted  that  it  exhibits  the  same  essential 

•ff 

character  as  the  cylinder  t distribution.  Note  in  passing  that  the  t dis- 
tribution of  Figure  19  is  locally  disturbed  by  the  presence  of  the  gill  hole, 
as  anticipated  in  paragraph  3,  page  21. 

This  reference  to  the  cylinder  problem  clarifies  somewhat  the  difficulty 
with  the  top  surface  boundary  condition  in  the  research  problem,  or,  if  pre- 
ferred, the  difficulty  with  the  contribution  from  component  problem  (iii). 

If  we  imagine  the  blade  to  be  (crudely  but  essentially)  modelled  by  a hollow 
unconstrained  cylinder  of  finite  length,  the  top  surface  boundary  condition 
used  in  (iii)  represents  condi tiors  near  the  midplane  of  a long  cylinder. 
Nearer  the  ends,  or  if  the  cylinder  is  shorter,  e.g.  only  two  ot  three  outer 


diameters  long,  the  top  surface  condition  up  = 0 leads  to  values  of  axial 


tension  on  the  inner  surface  and  compression  on  the  outer  surface  which  are 
higher  than  probably  realized.  It  is  therefore  difficult  to  judge  the  validity 
of  the  contribution  of  component  problem  (iii)  for  a specific  gill  hole,  de- 
pending on  its  location  relative  to  the  ends,  i.e.  top  end  and  throat,  of  the 
blade.  Nevertheless,  the  nature  of  the  contribution  of  this  component  to  the 
stress  at  point  A is  similar  to  the  contribution  of  the  centrifugal  field, 
whereas  the  contributions  of  component  problems  (i)  and  (ii)  were  comparatively 
little  understood. 
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6.  SOME  FINAL  REMARKS 

We  think  it  is  appropriate  to  close  this  investigation  by  briefly  stating 
what,  in  our  view,  are  the  more  important  contributions  arising  from  it  and  by 
summarizing  what  we  think  to  be  of  questionable  validity  and  why.  In  addition, 
a few  remarks  will  be  made  regarding  suggestions  for  improving  reported  data 
through  further  research. 

Within  the  continuum  theory  of  linear,  homogeneous,  isotropic,  steady- 
state  (classical)  thermoelasticity,  this  research  effort  has  produced  four 
separate  computer  programs  capable  of  performing  reliable  numerical  analyses 
for  a wide  variety  of  plane  and  three-dimensional  problems  governed  by  the 
theory.  This  computing  capability  is  quite  contemporary,  competitive,  and  in 
some  cases  better,  we  judge,  than  alternative  numerical  methods  for  problems 
of  thermoelasticity.  Indeed,  the  amount  of  interest  in  the  BIE  method  in  gen- 
eral, and  our  specific  programs  in  particular,  shown  by  the  academic  and  in- 
dustrial communities  here  and  abroad  is,  perhaps,  witness  to  that  fact.  The 
existence,  therefore,  of  the  four  programs,  motivated  by  and  created  largely 
for  the  solution  of  the  primary  research  problem,  is  probably  the  most  import- 
ant outcome  of  this  investigation. 

We  should  emphasize  that  these  programs  may  be  used  effectively  for  a 
broad  class  of  problems  in  elasticity,  with  or  without  thermal  effects  or  body 
force.  Clearly  the  data  and  problem  descriptions  for  the  present  work  were 
quite  specific.  However,  the  programs  may  be  used  for  analysis  using  data  for 
other  problem  descriptions  and  thermomechanical  configurations,  as  needed,  in 
whatever  discipline  classical  heat  conduction  and/or  (thermo)  elasticity  prob- 
lems occur. 


Regarding  our  "solution"  of  the  primary  research  problem  and  the  described 
related  problems,  we  feel  that  there  is  sufficient  uncertainty  that  the  turbine 
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blade  material  is  modelled  accurately  by  classical  thermoelasticity  to  relate 
all  of  our  results  to  "reality"  only  with  prudent  reservation.  That  point 
being  taken,  it  is  important  to  recognize  that  output  data  is  no  more  reliable 
than  input  data  and  herein  is  the  source  of  some  uncertainty  regarding  some 
of  our  results.  Specifically,  we  were  prepared  to  enter  thermal  data  (e.g. 
in  Eq.  (17))  with  considerable  detail,  but  the  availability  of  such  data  from 
the  blade  manufacturer  was  fairly  crude.  This  was  less  true  regarding  the  x,y 
variation  of  such  data  as  input  into  our  plane  analyses  but  sufficiently  true 
for  our  three-dimensional,  non-isothermal  analyses  to  cause  concern  about  the 
meaning  of  the  output  data  beyond  its  being  representative,  certainly  in  the 
"ball  park",  and  reasonably  self  consistent  whenever  self  consistency  checks 
were  available  and  made.  On  the  other  hand,  our  isothermal  solutions,  e.g. 
the  centrifugal  field  results  (Figures  15  and  16),  we  judge  to  be  excellent, 
granting,  of  course,  whatever  uncertainties  accrue  to  the  two  input  load 
models  (boundary  conditions)  used.  Further,  component  problem  (i),  leading 
to  the  results  of  Figures  13  and  14,  are,  we  judge,  similarly  good  despite 
the  dependence  of  the  shown  results  on  thermal  data.  This  is  true  since  the 
thermal  data  at  issue  for  problem  (i)  is  exclusively  plane  and  given  by  the 
manufacturer  in  probably  sufficient  detail  for  the  purpose  at  hand.  Added  to 
the  thermal  data  uncertainty  for  component  problems  (ii)  and  (iii),  and  hence 
for  the  full  thermoelastic  problem  (the  "sum"  of  problems  (i),  (ii),  and  (iii)), 
is  the  uncertainty  previously  discussed  at  length  regarding  the  horizontal 
surface  (Figure  10)  boundary  conditions. 

Despite  the  above-mentioned  difficulties  with  certain  of  the  data,  our 
analysis,  broken  down  as  described,  gave  us,  and  hopefully  will  give  the  read- 
er, valuable  insight  into  the  relative  roles  played  by  the  various  components 
of  the  thermal  field  and  the  centrifugal  field  in  producing  stress  at  the 


41 


H-T  

In 

I • i 

feather  edge  of  a gill  hole. 

By  way  of  further  research,  a renewed  attack  on  the  primary  research  pro- 
blem model  could  be  profitably  undertaken  with  the  availability  of  more  de- 
tailed input  thermal  data.  To  obtain  such  data  for  a more  detailed  descrip- 
tion of  the  parameters  hc  and  6Q  would  require  fluid-flow/heat- transfer  re- 
• search  beyond  the  scope  of  the  present  investigation.  A perhaps  better  set 

of  mechanical  boundary  conditions  than  those  used  for  Figure  10  could  be 

I ' 

obtained  from  a three-dimensional  thermoelastic  analysis  of  the  entire  blade 
minus,  of  course,  the  individual  gill  holes.  Here  again,  however,  the  axial 
variation  of  temperature  would  have  to  be  provided  from  the  manufacturer.  A 
fixed  throat  portion  of  the  blade  with  the  other  end  free  would  perhaps  pro- 
vide realistic  mechanical  boundary  conditions  for  such  a model.  From  a sol- 
ution for  this  overall  model,  it  would  be  possible  to  determine  a possibly 
better  "top-surface"  boundary  condition  for  the  model  of  Figure  10  plus  bet- 
ter tractions  on  the  vertical  surfaces  there  with  some  axial  (z)  variation. 

Finally,  the  broader  questions  concerning  the  thermoelastic  model  it- 
self, e.g.  the  questions  of  inelastic  effects,  low-cycle  fatigue,  crack  pro- 
pagation, temperature- dependent  material  properties,  are  all  relevant  to  the 
described  physical  turbine  blade  configuration  and  should  be  addressed  as 
feasible  in  future  investigations. 
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Thermal  Stress  Distribution  on  Piece  of  Turbine  Blade 


Figure  9.  t Distribution  over  Blade  Cross-Section 
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Distribution  for  Full  Thermoelastic  Problem 
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Figure  13.  Stress  Distribution  for  Component  Problem  (i) 


Figure  14.  Hoop  Stress  Pattern  for  Problem  (i) 


15.  Stress  Distribution  for  Centrifugal  Input  a 


Figure  16.  Stress  Distribution  for  Centrifugal  Input 
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